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Abstract 

We consider in this paper two methods to differentiate functionals of the kind Jg n da when 
/ depends only on the geometry of d£l. The first method is based on shape optimization with 
a Reynolds formula for surfaces. The second consists in a level set approximation on a volumic 
domain. The main contribution of this paper is to prove that these two approaches lead to the 
same geometric result when / depends on the normal and the mean curvature and allow to 
consider general functional depending on the Gaussian curvature. 

Introduction 

An energy functional J depending on a domain $7 is associated in a large range of problems 

J (SI) = [ /(O) da 
Jon 

In order to minimize this energy with a gradient method we have to differentiate J(O) with 
respect to the domain. Shape optimization is now classical. It was introduced originally in appli- 
cations where / is solution of a PDE [?, ?, ?, ?]. We focus here on geometrical shape optimization, 
in this case / depends only on the geometry of £1. These energies arise in a lot of applications, 
for example anisotropic mean curvature flow [] , image processing [] and shape equilibrium of 
vesicles []. 

A way to perform the geometric shape optimization is to introduce a parametrization of the 
surface and to compute the derivative in the normal direction. This approach has been used 
by Willmore for the mean curvature [?] and by Steigmann for the Gaussian curvature [?] and 
leads to a geometrical result. However complex tools of differential geometry as the covariant 
derivative are needed in the computations. To get round of these difficulties we propose two 
approaches based on a level set representation of the functional. 



The scope of the paper is the following. In the first section by introducing shape operators and 
curvature. Then we prove some important lemmas and integration by parts. The second section 
is devoted to shape optimization of geometrical functionals with a Reynolds formula for surface. 
The differentiation of functional depending on the normal and the mean curvature are carried 
out and the result is given in a geometrical form. The third section deals with another approach 
based on a volumic approximation with level set functions. Doing so we rewrite the derivative for 
functionals depending on the normal, mean curvature and Gaussian curvature (with a distance 
function). We find the results of the second section. 



1 Notations and preliminaries 

In the sequel, O will be an regular open connected subset of M 3 . Note that its boundary dfl is 
a smooth closed surface of M 3 . We denote by ■ the euclidean scalar product on M 3 and | • | the 



associated norm. The vectors and matrix will be written in the canonical basis of R 3 . For two 
vectors a and b, a®b stands for the tensor product that is [a(8>6]y = a>ibj. For two matrix A, B we 
note A: B = J\j MjBij their scalar product. We note Tr(A), Tr(Cof(A)) = ±(Ti(A) 2 -Tr(A 2 )) 
and det(A) the three invariants of a matrix A in M3(R). We denote by V, div and A the usual 
gradient, divergence and laplacian, respectively. 

We denote by re a unitary smooth extension on R 3 of the normal to the surface dfl. Taking the 
gradient of |n| 2 = 1 gives 

[Vnfn = 0, ([Vn] re) ■ re = 0. (1) 

These equations will be used extensively in this paper. Let us first introduce some basics tools 
of differential geometry as shape operators and curvature. Then we prove some lemmas and 
formulas to integrate by parts. We conclude this section by introducing the framework of level 
set methods. 

1.1 Shape operators 

The shape operators generalize classic differential operators on a surface. The reader may consult 
[?, ?] for a more detailed exposition of this material. The tangential gradient of a smooth function 
/ : R 3 — ► R is defined by 

V 9 n/ = (/ - n ® n)V/ = V/ - (V/ • n)n. 

The tangential gradient only depends on the values of / on the surface d£l. A simple consequence 
of the definition is that Vgn/ • re = 0. The tangential gradient of a smooth vector field v : R 3 — > 
R 3 is the tangential gradient of each component 

[V 9 nv] = [Vv](I -n®n). 
The tangential divergence of the vector field v is the trace of its tangential gradient 

divdn(v) = TrQVanf]) = div(v) - ([Vv] re) • re. 
The Laplace-Beltrami operator on d£l of / is the tangential divergence of the tangential gradient 

&dnf = div 9n (V 9 n/)- 

These operators are defined on R 3 but they have a geometric interpretation on the surface dQ. 

1.2 Curvature 

The curvature is defined as the variation of the normal in the tangent plane and we denote by 
n\ and K2 the eigenvalues of the associate endomorphism dn. By definition, the mean curvature 
H is Ki + K2 and the Gaussian curvature G is k\K2- The extension of dn on the whole space is 
the shape operator associate to the normal [Vgn n ]- The principal curvatures k\ and K2 are then 
the nonzero eigenvalues of [Van"-]- Therefore 

H = Tr([V 9 nn]) = div 9 n(n), G = Tr(Cof ([V an n])). (2) 

If the extension of n is unitary we have 

H = TrQVre]) = div(re), G = Tr(Cof([Vre])). (3) 

The proof is based on (pQ). Throughout this paper, we assume that the extension of the normal 
is unitary. 



1.3 Lemmas 

The Leibniz rule to derive the product of functions may be generalized shape operators. Consider 
two smooth functions /, g : K 3 — ► R. Following ([2]) we get 

divan (fn) = f div dn (n) + V dn f • n = fH. (4) 
The operator P = I — n<gm is a symmetric projector so 

Van/ • V 9 = V/ • Van 5 = V dn f • Vans- (5) 
The normal is unitary so n = A for some vector u. We have the following identity 

[Vn] n = A (V\v\ - ( VH • ^] = (6) 
\v\ V V M / M / \v\ 



We now prove some useful lemmas for the next sections of the paper. 
Lemma 1 

div([Vn] n) = VH ■ n + H 2 - 2G. 

Proof : The property ([3]) implies 

divflVn] n) = VH ■ n + Tr([Vn] 2 ) 
The lemma is proved with the identity Tr([Vn] 2 ) = k 2 + k| = H 2 — 2G. 



Lemma 2 

(V(V an /)n)-n = -V/-([Vn] n). 



Proof : Taking the gradient of the equality Van/ ■ n = and multiplying by n gives 

([V(Van/)] T n) • n + ([Vn] T Van/) • n = 0. 
The lemma is proved combining ([6]) and ©. 

Lemma 3 

(V ([Vn] n)n)-n = -([Vn] n) • ([Vn] n). 

Proof : Taking the gradient of (pQ) gives 

[V([Vn] n)] T n + [Vn] T [Vn] ra = 0. 
Multiplying by n the result prove the lemma. 



1.4 Integration by parts 

Let be Q an open set of R 3 . Consider two smooth functions /, g : M 3 — > R and a smooth vector 
field f : M 3 — > M 3 . Assume that all these fields vanish on dQ. We have the following integration 
by parts for shape operators 



/ divan (u) dx = - Van/ • v dx + / div(n (gin) ■ v fdx (7) 
Jq Jq 

The proof is a simple consequence of the Stokes formula and the definition of shape operators. 
We can also integrate by parts shape operators on a closed surface <9$7 

f Vaaf-v = -f f div dn (v) da + [ Hfvnda (8) 
Jan Jan Jan 

The proof of this result can be found in [?] p 55. We provide another demonstration with a 
volumic regularization in the third section. Compared to the classical Stokes formula in M 3 for 
usual differential operators, an extra term involving the mean curvature of the surface appear. 
A consequence of flSJ) is the following 

/ f^an9 da = gA m f da (9) 
Jan Jan 

1.5 Level set formulation 

Consider a smooth function 4> : R 3 — ► M whose level set represent the smooth surface dQ 

on = {x e R 3 I 4>{x) = 0} 

Assume that |V</>| > in a neighbourhood of {(ft = 0}. Consequently, the unitary normal is 
n = |y^i- If ^ is a distance function ie \V(fi\ = 1 then [Vn] is symmetric and ([1]) becomes 

[Vn] n = (10) 



In the framework of level set methods, the geometrical functional /(^) da may be written 

J(0) = / f[<f>] da (11) 
J{<f>=0} 

where the notation f[(f>] stands for a function / depending on eft and its derivatives. A natural 
idea to perform the shape optimization of J(4>) is to differentiate with respect to </>. However, we 
prefer to use an equivalent method which consist in transporting eft (which from now on depends 
on time) with a velocity field u(x, t) 

(f) t + u-V4> = (12) 

Therefore, the shape optimization of J (eft) is performed by differentiating with respect to the 
parameter t. To get round of the problem of differentiating on a domain depending of 4> we 
propose two methods. The first one consists in applying a Reynolds formula for surfaces and the 
other one is based on a volumic regularization of the functional. 



A 



2 Optimization with a Reynolds formula for surfaces 

A proof of the Reynolds formula for surfaces is carried out in the annex. Integrating by parts 
(1301) witn © gives 



{J{4>))t = (M)t + div(/[0]n) u ■ n da. (13) 

J{4>=0} 

We get the formula obtained in the framework of shape optimization develop by Murat and 
Simon [?]. We focus now on the particular case where / depends on the normal and the mean 
curvature. 

2.1 Functional depending on the normal 

We first compute the derivative of the normal with (fT2l) and ([6]) 



{nmt = w\~ ww m = — m — = m(u ' n) " ([ ] n)u ' n - (14) 

We get the result obtained by shape optimization in [?] p 70. 

The normal is unitary so it is natural to introduce a smooth function / : R 3 — ► R homogeneous 
of degre one and the functional 



J n (cP) = / f(n(</>)) da. 
J{4>=0} 

By ()13p and the derivative of the normal (|14p we get 



(J n ((f>))t = / div(/(n)n) u • n + V/(n) • (-V«i(« • n) - ([Vn] n) u ■ n) da. 
J{4>=0} 

Using the symmetry relation © and integrating by parts with flSJ) 

-/ V f(n)-V dn (u-n)da = - / V dn f{n)-V 9n (u-n)da= div m (Van f(n))u-n da. 

J{4>=0} J{4>=0} J{<i>=Q} 

By © 

div(/(n)n) = /(n) div(n) + V(/(n)) • n = f(n)H + ([Vn] T V/(n)) • n. 

Therefore 



(J n (</>))t = / (/(n)fr + diven(Van/(n)))«-nd<7. (15) 

J{<l>=0} 

The relation f(tn) = tf(n) imply V/(n) • n = f(n). Hence by 

divan (Van /(n)) = divan (V/(n)) - divan((V/(n) ■n)n) = div S n(V/(n)) - f(n)H. 
Finally 



(J n (</>))t = / divan(V/(n))n-nda. (16) 

We get the anisotropic curvature Hf = divan {nf] with the anisotropic normal rif = V/(n) as 
in [?] p 192. 



2.2 Functional depending on the mean curvature 

Consider a smooth function A : R — ► R and the functional 

J H {4>) = [ A(H(4>)) da. 
J{4>=0} 

The integration by parts formula ([8]) involve shape operators so we choose the definition ([2]) of 
the mean curvature. Following fjl3[) we get 

(J H (<P)) t = [ dw(A(H)n)u-n + A'(H)(H(ct>)) t da. 

The derivative of the normal (|14p gives 

(H{<j>)) t = div(-V d n(u ■ n) - ([Vn] n) u ■ n) - (V(-Van(« • n) - ([Vn] n) u ■ n) n) ■ n 

— ([Vn] (— Vqq(u ■ n) — ([Vn] n) u • n)) • n — ([Vn] n) • (—Vgn(u ■ n) — ([Vn] n) n • 

We put the two first terms in a tangential divergence and the third one vanish by (P) . Developping 
the divergence gives 

(H((f)))t = -Asn(u ■ n) - div9o([Vn] n) u ■ n + ([Vn] n) ■ ([Vn] n) u ■ n. 
Following lemma Q] and lemma [3] 

div 9n ([Vn] n) = div([Vn] n) - (V([Vn] n)n) ■ n = VH -n + H 2 -2G+ ([Vn] n) • ([Vn] n). 

We have also div(A(iJ)n) = A(H)H + A'(H)X7H ■ n. After integrating by parts the laplacian 
with we finally get 

(J H (<P))t= [ (A(H)H-A dn (A\H))-A\H)(H 2 -2G))u-nda. (17) 
J{4>=0} 

We get the result of Willmore [?] which is geometrical because it involves only the curvature and 
shape operators. For a more general functional depending on the Gaussian curvature, complex 
computations gives the same result as the one developed in next section devoted to optimization 
by volumic approximation. 

3 Optimization with a volumic approximation 

The following lemma is the key of method of optimization by volumic approximation. 

Lemma 4 Let be £ a continuous positive cut-off function with support in [—1, 1] and J" R £ = 1. 
Assume that f,cp are regular and |V</>| > 0, then 

U4>) := [ fl<P}\V<P\-( (7) dx — / M da 



Proof : The co-area formula is 

I f(x)dx= [ [ /(aOlV^r 1 dadv (18) 

J {\4>(x)\<e} J-e J{<t>{x)=v} 

This formula is quite natural since the volume {\<f>\ < e} is computed integrating the area of 
{4> = v} with a factor |V</>| -1 which correspond to the spacing between the level set of (p. We 
refer to [?] p 118 for a rigorous proof. For a continuous function g and under the hypothesis 
on C we have ff (7) <?( r ) dr — > g(0). Applying the co-area formula (fT8|) to the continuous 

function g(r) = JV^* , /(x) do" prove the lemma. 



Remark 

This lemma allows us to get the integration by parts ([5]). Indeed, the integration by parts for 
shape operators ((TJ), the property Vgafi = and ([6]) give 

J^div dQ (v) |V0|^C(~) dx = J Hvn\V<t>\~ £ Q(^j dx. 
Taking e — ► with the lemma [Hand replace v by fv gives the result. 

3.1 General functionnal 

We compute the derivative of J e {4>) with respect to t 

(Je(4>))t = J(M) t mk (^j dx+J /[0]|c (£) j|||-v& dx+J mm^c (P) & dx. 

Integrating by parts the second term gives 



Jq ^ (7) ^ - = - j Q *° 7< (!) * (7) 

Hence (112p implies 

(JM)t = J ( (f[<t>]) t + div (/[0]n) « • n ) |Wl§C ( dx. (19) 

We get the Reynolds formula for surfaces (|13|) taking e — ► with the lemma 01 We focus now 
on the case where / depends on the normal and the curvature. 

3.2 Functional depending of the normal 

Consider a smooth function / : R 3 — > M homegenous of degree one and the functional 

Je(<P)= f /(n(0))|V^c(j) dx. 



Following (|19[) and the derivative of the normal (jl4[) we obtain 

(TO))t = y div(/(n)n)«-n|V^|^c(|) + / V/(n) • Van&^C (7) 
Using the symetric relation ([5]) and integrating by parts gives 



J V/(n)-V 9n ^~C (jj dx = J Van/(n)-V^ic (7) dx = ~ J Q div (Van/(n))^C (j) & 
(fT2l) implies 

(J2{4>)) t = J div(/(n)n + V 9n /(n))u-n|V0|^C^J dx. (20) 
This formula may be written in another form. Indeed with lemma [2] 



div(Van/(n)) = div an (V^/(n)) + (V(V 9Q /(n)) n) • n = div sn (V 9 n/(n)) - V/(n) • ([Vn] n). 

We have also div(/(n)n) = f(n)H+ ([Vn] T V/(n)) -n. Therefore we get the previous result (|15[) 
taking e — ► with lemma [H 

3.3 Functional depending of the curvature 

Consider a smooth function g : M3QR) — ► K and the general functionnal 

J e (^) = y 5(V(n(0)))[V^c(j) dx. 

We denote by [Vg] the gradient of 5 at the point [Vn]. Following (|19p and the derivative of the 
normal £0 



(J E (0)) t = jTdiv(^([Vn])n)«-n|V0|ic(j) dx + ^[V<?] : V (^f^) |V<^C (j) dx. 
Integrating by parts the second term gives 



The integration by parts ([7]) for shape operators gives 



dx. 



mm 



j div(<?([Vn])n) u ■ n |V<^[— C, ^— J dx + J div ^[Vg]|V(/>[ — £ ^— • div(n ® n) u • n fix 
- J Q div 9 n (^ div ([ V ^l V ^l^c(^))) u-n|V</>|dx 
Developping the divergence 

div ([vg]m-c (*)) = div([v 5 ])|v^i-c (£)+\yg] viv^i-c (£)+\yg] m\c 



£ ' \ £ J J £ \£ J ' ' ' ' £ ' \ £ J " ' ' - 1 \ - 

The relation V 9 n(</>) = implies divg^ (&j = div 9 n(^)^C J f° r an vector v. We obtain 
finally 



(JM)t = I div( 5 ([Vn])n) u ■ n |V0|-C ( -) dx - [ div 9f2 (div([Vg])) u ■ n \Vcj>\-( I - ) dx 



Q c \ t / JQ 



£ \£ I .In £ \ £ 



+ ^div([V 5 ])-div(n®n) M .n|V0|ic(^) dx + / g (Vs] • 
+ / ([Vg] V</)) • div(n® n)^C' ( - ) w • n |V<£| dx. 



div(n ®n)u-n IV^I— £ 

£ 



u ■ n |V^| 



In order the simplify this complex expression in the case of a functional depending on the mean 
and Gaussian curvature we consider 

g{w) = F(Tr(uO,Tr(CofH)). 

where F:lxl — ► R is a smooth function. We denote Fh and Fq the derivative of F with 
respect to his first and second variable, respectively. We have 

[Vg}([Vn]) = al + PlVnf. 
where a = Fh + HFq and = —Fq- 

The two last terms involving in the expression of (J £ ((f)))t seem annoying. However using (JTJ) 
and (|U) these terms compensate. 

We now simplify the fourth and fifth terms involving Hy-ffi ■ Following we have = 
[Vn] n + V |^|'" n. Hence with ® and © 

- divan ((a/ + /3[Vn] T )^^ = - div 9n ((a/ + /?[Vn] T )[Vn] n) - aff V| ^J' W 
We have also 

iTn V|V</)|\ _ wr _ , , V|V0|-n 



(a/ + /3[Vn] J ) , ' ,, -div(n®n) = (a+/?if) ( [Vn] n)-([Vn] n)+aH ' ,', — +/3([Vn] [Vn] n)-([Vn] n). 
|V0| / 1 |V0| 

We denote by I £ the expression of (J £ (4>))t without the integrals and the terms u ■ n|V0|^C 
Therefore 



I £ = div( 5 ([Vra])n) - div 9n (div(a/ + 0[Vnf)) - div a n ((a J + 0[Vnf) [Vn] n) 

+ div((al + (3[Vn] T )) ■ div(n n) + (a + /3if)([Vn] n) • ([Vn] n) + /?([Vn] [Vn] n) • ([Vn] n) 

We propose to compute the contributions of al and /?[Vn] T in order to obtain the derivative of 
a functional depending of the mean and the Gaussian curvature. 

Functional depending on the mean curvature 

We note R a the expression J e with = and without the term div(<?( [Vn])n) 

R a = - divan(Va) - divan (a[Vn] n) + Va • (Hn + [Vn] n) + a([Vn] n) • ([Vn] n). 



a 



(|3J) implies divaQ(Va) = Ag^a + HVa ■ n. Following lemma El ([5]) and ([6]) 

— divgn(a[Vn] n) = — adiv([Vn] n) — a([Vn] n) • ([Vn] n) — Va • ([Vn] n). 
We finally get 

i? Q = -A dn a - a div([Vn] n). (21) 

Consider the functional 

j? (0) = J ^(iJ(^))|V0||c (£) ^ 

which correspond to a = ^4'(i?) and /? = 0. The relation d\w(A(H)n) = A(H)H + A'(H)VH ■ n 
and the lemma [T] gives 

(Jf = | ( A(fT)lT - A 9n (^(#)) - A'(tf )(if 2 - 2G) ) « • n |Wl|c (~) 
We get the result (|17|) taking g — ► with the lemma HJ 

Functional depending on the Gaussian curvature with a distance function 

We choose a distance function for <p to simplify I £ . We note the expression of I £ with a = 
and without the term div(g([Vn])ra) in the case of a distance function. Therefore, by (|1U|) 

RP = -div^(div(/?[Vn] T )) + div((/3[Vn] T )) • ffn. 

Following ([TO]) we get div(/?[Vn] T ) -f^n = f3HVH-n. We have also with (JT]), ([3]) and the property 
of symmetry ([5]) 

- div 9Q (div(/3[Vn] T )) = -pA dn H - (3H(VH ■ n) - V dQ H ■ V an /3 - div^([Vn] T V dQ (3). 
The relation divgfi(Ab) = (1wqq(A t ) ■ b + A T : V ' g^b and the symmetry of [Vn] gives 

-div m ([Vn] T V 9Q /3) = -div m ([Vn] T ) • V dn (3 - [Vn] : Vg n {V dn (3). 

Taking the gradient of (pQ) and multiply by n gives (V([Vn] T )n) • n = -[Vn] T [Vn] n = 0. 
Therefore divg^([Vn] T ) = VH. The symetric relation © and the relation [Vn] = [V^n] valid 
for a distance function gives 

R P d = -[3A dn H - 2Vg n H ■ V dn [3 - [X7 dn n] : V 9n (V sn /3). (22) 

Consider the general functional 

J?> G (0) = J Q F(H(0),G(<p))mk (f) dx 

Let J d the sum of R a $Z$, R% J22J and the term dw(F(H,G)n) with a = F H + HF G and 
/3 = —Fq. In case of a distance function we have R a = —Ag^a so 

J d = dw(F(H, G)n) - A dn (F H + HF G ) + F G A dn H + 2Vg Q H ■ V dn F G + [V an n] : V dn {V dn F G ). 



1 n 



Developing the surfacic laplacian 



-A dn (HF G ) + F G A dn H + 2V dQ H • V dn F G = -HA dn F G . 
Using the lemma Q] with a distance function 

div(F(F, G)n) = F(H, G)H - F H {H 2 - 2G) + F G VG • n. 
Taking e — ► with the lemma H] gives 

{J H > G {4>)) t = [ (F(H,G)H -F H {H 2 -2G)-A dn F H )u-nda (23) 

J{4>=0} 

+ f ( F G VG • n — HA dn F G + [V dn n] : V dn (V dn F G ) )u-nda (24) 
J{4>=0} 

Assume the surface do not change of topology. The Gauss-Bonnet theorem claims that /r^_ i G da- 
is constant. Applying the formula (|24j) with F(H,G) = G gives (Fjj = and F G = 1) 
I{4,=q}{GH + VG ■ n) u ■ n da = 0. This relation is valid for all u so VG • n = —GH on 
{0 = 0} 

We obtain finally 

(J H ' G ( ( t>)) t = [ (F{H,G)H -F H (H 2 -2G)-A dn F H )u-nda (25) 

J{4>=0} 

+ f ( —F G GH — HA dn F G + [V dn n] : V dn (V dn F G ) )u-nda (26) 
J{4>=0} 

We get the result of Steigmann [?] which is geometric because it involve only shape operators 
and curvature. 

4 Conclusion 

In this paper two methods of shape optimization had been investigated for functionals depending 
on the normal and the curvature. The first one is based on a Reynolds formula for surfaces and 
integration by parts forces us to introduce shape operators. The second one, based on volumic 
approximation with level set leads to simple computations but some efforts are needed to put the 
results on a geometrical form. These two approaches give the same result which is geometrical 
for functional depending on the normal and the mean curvature. The second methods allow us 
to consider more general functionals depending on the Gaussian curvature. These formulas will 
be useful to study the influence of the Gaussian curvature in the shape equilibrium of vesicles. 

5 Annex 

To be self contained we present here a demonstration of the Reynolds formula for surfaces based 
on change of variables on surfaces. 



1 1 



5.1 Change of variables in surface integral 

Let be U an open set of R 2 and a smooth function $ : R 3 x R — ► R. Consider X : U — ► R 3 a 
smooth parametrisation of the surface OQq and Y : U X R — > R 3 a parametrisation associate 
to the surface <9Qj = $(9f2o>i)- Consider also a smooth function / : R 3 x R — > R. We have 

f(x,t)da= [ f(Y(u 1 ,u 2 ,t))\Y 1 AY )2 \du 1 du 2 . 
The relation Yj = [V3>] Xj gives for all vector v 

(Yi A Y 2 ) • ([V$] u) = det([V$] [W] X 2 , [W] u) = det([V$])det(Xi,X 2 ,u). 
Therefore 



FiAF 2 = det([V$])[V$n T Xi A X 2 . (27) 



The relations uq = H^rVn and n = ly 1 ^ 2 implies |Yi A Y 2 | = | Cof([V<3?])no||X i A X 2 | so 



f(x,t)da = /(*(o,t),t)|Cof([V*](a,*))no(a)|d7o- (28) 

<9f2i J9Q 



We have also from (12 



_ Cof([V^])n 

|Cof([V$])no|' 1 j 

5.2 Reynolds formula for surfaces 

Consider a surface <9f2t evolving in a smooth velocity field u : R 3 x R — > R 3 and / : R 3 x R — ► R 
a smooth function. Then we have the following Reynolds formula for surfaces 

1([ f(x,t)da)=[ f t + u-Vf + fdiv dn (u)da. (30) 
dt \JdQt J JdQt 

Proof : The main idea consist to write the integral on a fixed domain d£lo and differentiate 

with respect to t. Consider the caracteristics 3> associates to the velocity field u ie ^|(a,t) = 

u(&(a,t),t) and denote dQt = $(dQ,o,t). The change of variables x = $(a,t) ([28j) gives 

f(x,t) da= [ f(^a,t),t)det([V$(a,t)})\lV<S>(a,t)}- T no(a)\ da . 
idsh Jan 

Differentiating with respect to t gives three terms. The caracteristics equation gives for the first 
term (f(<&,t))t = ft + u- V/. The caracteristics equation gives also 

[V$] t = [Vu][V$] (31) 

The derivative of the determinant (det(A(£)))' = det(A(t)) Tr([A(t)] _1 A'(t)) and ([3TJ gives for 
the second term 

(det(V*)) t = det(V$)div(u). 

For the third term, the derivative of the inverse (A' 1 )'^) = -A^^A'^A' 1 ^), flH} and (f29l) 
gives 

|[V$]- T n | t = (([V^]- T ) t n ) ■ j^j^ n ° = - (N n) • n |[V#]~ T n |. 

|[V5>J x no| 

The change of variables = $(a, t) and the definition of tangential divergence prove the lemma. 
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